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Abstract. We characterize the flickering observed in the optical lightcurve of the Intermediate Polar system V709 Cas by 
determining its position in the a-S as in the Fritz and Bruch (1998) classification scheme. S represents the strength of flickering 
at a given timescale, while a describes the energy distribution of the flickering at different time scales. Here a is independently 
obtained with both the wavelets and the Hurst R/S analysis. The flickering shows self-similarity in the time scale ranging 
from tens of minutes down to 10 seconds with stochastic persistent memory in time, a and S appear anticorrelated. In the a-Y, 
diagram V709 Cas falls in the region of magnetic systems. Since V709 Cas shows the spin period of the magnetic WD only in 
the X-ray but not in the optical, we conclude that this method can be used to characterize CV subtypes especially when their 
classification is uncertain. 
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1. Introduction 

Cataclysmic Variables (CVs) are binary systems in which a 
late-type secondary star fills its Roche lobe and transfers mat- 
ter onto a white dwarf (WD) primary. The WD is surrounded 
by an accretion disk, unless its magnetic field is strong enough 
to partially or totally control the accretion geometry. 

CVs are classified into three major groups: nova-like (NL) 
systems, classical and recurrent novae (CN, RN), and dwarf 
novae (DN). In particular, nova systems must host CNO WDs 
with masses larger than ~ O.6M so that the accreted material 
can be cyclically ejected through nova outbursts (Livio 1992). 
CVs can also be classified according to the strength of the mag- 
netic field of the primary. Thus, we have the non-magnetic 
systems, the Intermediate Polar systems or DQ Her systems 
(hereafter IPs), and the Polars or AM Her systems. In the non- 
magnetic systems accretion onto the WD occurs from the last 
stable orbit of the inner disk. In IP systems the magnetic field 
disrupts the inner regions of the accretion disk within the so- 
called Alfven radius and accretion mainly occurs via accretion 
curtains onto the magnetic poles of the WD. In Polar systems 
the disk is totally disrupted by the strong magnetic field of 
the WD and the accretion stream is directly conveyed on the 
WD magnetic poles. In both Polars and IPs the infalling matter 
forms a strong shock above the magnetic poles of the compact 
star that mainly radiate in the X-ray domain and eventually 
producing cyclotron radiation in the optical/IR. 

At any given orbital period, the luminosities of NL systems 
and of quiescent novae are systematically brighter than those 



of DN, suggesting that the former ones are powered by larger 
mass transfer rates (Mt) from the secondary and/or have hotter 
WDs preventing the onset of the cyclical disk instability phe- 
nomena observed in DN. An exception may be represented by 
those CN and NL that are also magnetic Polars in which the 
strong magnetic field of the WD reduces the efficiency of the 
magnetic braking mechanism of the binary system and, con- 
sequently, the mass transfer rate and the accretion luminosity 
dWarner, 1995) . 

The lightcurves of Cvs may show a variety of periodic, 
quasi-periodic and/or erratic modulations. Flickering consists 
of short term variations on time scales from a few seconds 
to a few tens of minutes, with magnitude variations from a 
few tenths of magnitude to the magnitude scale. Flickering 
appears as a sequence of overlapping flares and bursts with a 
random variability in time, sometimes showing self-similarity 
at different time scales, i.e. the stochastic fluctuations in 
the detrended lightcurve are similar independently from the 
chosen time binning. Warner & Nather (1971) first tried to 
associate the flickering of U Gem and DN to the hot-spot 
in the outer rim of the AD although this ansatz was not 
confirmed in the observations of some other CVs like VW 
Hyi dWarner, 19751 |van A meronge net al., 1987| l, Z Cha, 
V2051 Oph, ( |Warner & O'Donoghue, 1987] >. The observations 
of AE Aqr and YZ Cnc (Elsworth & James 1982, 1986; 
James, 1987), instead, suggested that the origin of flickering 
might lie in the innermost regions of the AD or close to the 
surface of the accreting WD. Some models suggest that the 
instabilities in the mass accretion onto the WD (Bruch 1992, 
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1992a, 1992b, 1996) are due to MHD plasma turbulence 
in the accretion flow. Recently, the flickering of CVs was 
studied in several systematic high speed photometry cam- 
paigns (IWoudt & Warner 200 1 1 IWoudt & Warner 20021 

IWoudt & Warner 2002bl IWoudt & Warner 20031 

IWoudt & Warner 2003bl IWoudt et al. 20041 IWoudl et al. 20051 
IWoudt & Warner 20081 

Flickering has been observed both in non magnetic and 
magnetic CVs. Magnetic CVs usually present flickering in the 
optical and sometimes in the X-ray region of the spectrum. 
One general property of flickering is the correlation between 
its amplitude and the luminosity of the "quiet primary". The 
"quiet primary" is defined as the sum of the luminosity of the 
WD and the luminosity of that part of the accretion process 
(disk or flow) which is not involved in the flickering activity 
(Zamanov & Bruch, 1998). Information can be also obtained by 
comparing the flare rate and the ratio of the maximum and the 
mean flux observed to that of the quiet primary. In some objects 
the flickering involves a large fraction of the total light emitted 
eventually dominating the optical light curve. The modalities 
of the flickering observed may vary from object to object but 
some common features have been observed in CVs of the same 
sub-type. Bruch et al. (Bruch 1992) proposed a classification 
method of the flickering showing that different subtypes of CVs 
tend to occupy specific regions of the a-1, plane, where a is the 
parameter related to the energy distribution of the flickering at 
different time scales and X represents the strength of the flick- 
ering for a given time scale. 

In this paper we derive the a parameter for the IP system 
V709 Cas both using the wavelets and the Hurst R/S analy- 
sis (Hurst et al. 1965) and discuss the two methods. While 
the wavelet analysis of flickering is dependent on the choice 
of the mother wavelet function, Hurst's R/S analysis bypasses 
this problem, giving an independent determination of a. In ad- 
dition, the value of the Hurst exponent gives the degree of per- 
sistence/antipersistence of the stochastic memory of the flick- 
ering, depicting the global behavior of the physics of the ac- 
cretion process. Differently from the X-ray band, this magnetic 
CV shows no spin modulation in the optical range and it is 
strongly flickering-dominated, without any precise and stable 
dominating frequency. V709 Cas hence represents an ideal test 
case for our study. 

In Sec. 2 we introduce the mathematical basis for both 
wavelets and Hurst's analysis and explain the link between 
these two methods. In Sec. 3 we characterize the flickering of 
this IP system by placing it in the a, S diagram. Sect. 4 draws 
the conclusions. 

2. Flickering properties and Hurst exponent 

In the past, many authors discussed flickering mainly in terms 
of statistical properties (Mumford, 1966, Robinson, 1973| 
IMoffet & Barnes, 1974[ |Zuckermann, 1971K Bruch (1992, 
1992b, 1996), Bruch & Griitter (1997). A further step was 
made by Fritz & Bruch (1998) who, following Scargle et al. 
(1993), introduced the use of wavelet transforms and of the 
scalegram. The scalegram is a plot of the logarithm of a ref- 
erence timescale t s vs the logarithm of the variance S(t s ), 



namely a measure of the variances of the wavelet coeffi- 
cients expressed as a function of the time scale. See also 
|Percival & Walden7"2 000 for more details. 

Flickering is a stochastic process in time whose properties 
can provide vital information about the driving mechanism be- 
hind the accretion process. One of the most important statisti- 
cal properties of flickering is the measure of the variance of the 
wavelet coefficients at different time scales 

k 

where c 2 , are the wavelet coefficients in which k is the time 

s,k 

index, N the number of measurements (equivalent to the total 
number of the wavelet coefficients) and s is the scale index, 
also called "octave", related to the time scale by t s = 2 s At. 
S'(t s ) is related with the total energy of the signal calculated 
at the time scale t s and summed over all the time bins, then 
normalized to N and to the integration interval At. The function 
S'(t s ) characterizes the statistical variations of the data record 
in the time interval [2 s Af, 2 I+1 At]. 

Fritz & Bruch (1998) and Zamanov & Bruch (1998) 
showed that the scalegrams of flickering are with a good 
approximation linear functions of t s for most of the CVs' 
lightcurves. This implies the presence of self-similarity with 
a power-law behavior S'(t s ) ~ s a within certain time scales. 
The stochastic fluctuations in the detrended lightcurve are sim- 
ilar independently from the time binning as described by l// 8 
noises, where / is the time frequency and /3 = a/2. 

The linearity of the scalegrams permits a simple 
parametrization of the flickering properties with two parame- 
ters, a and 2. a is the inclination of the scalegram with respect 
to the x axis, that can be then determined with a linear least 
squares fit to the scalegram points of the log S(t s ), log(f. s ) plot. 
a shows whether slow or rapid light fluctuations are dominating 
the stochastic time series and the flickering strength parameter 
I = logS(t s ) ■ t s , given by 

NAt 

S(t s ) = S'(t s )- - (2) 

*->s,k c s ,k 

is given by the measure of the variance of the wavelet coef- 
ficients for a given time scale. Some possible correlations be- 
tween the strength and the duration of the stochastic light vari- 
ations can be put in evidence by plotting a vs 2 in the so-called 
or-I! parameter space. 

An alternative robust statistical tool that can characterize 
independently from the wavelet approach the stochastic prop- 
erties of flickering (i.e. the parameter a) is Hurst's R/S anal- 
ysis. This method offers the advantage of being independent 
from any previous modeling of the problem, like the choice of 
a mother wavelet. The approach initially used by Hurst (1951) 
and Hurst et al. (1965) is based on the rescaled range analy- 
sis which consists in estimating the ratio between the range of 
the variations, R, and the standard deviation S derived from the 
analysis of all the sub-intervals of data for each equal parti- 
tioning of the full data record. Mandelbrot & van Ness (1968) 
and Mandelbrot & Wallis (1969) linked this method to a partic- 
ular class of self-similar random processes, named Fractional 
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Brownian Motions (FBMs). It was also shown that also short- 
run statistical dependence in pseudo-random sequences can be 
described with this technique (Gammel, 1997). Hurst found 
that natural phenomena follow the empirical power law 

R/S = (eN) H (3) 

where e is a constant, N the number of data present in each sub- 
interval, and H is the so-called Hurst's exponent. In particular, 
Hurst found that purely stochastic phenomena present a mean 
value H ~ 0.73 and a standard deviation S ~ 0.09. 

In dynamical systems, H characterizes the stochastic mem- 
ory in time of the process. Hurst exponents H > j indicate 
the persistence in time of a certain trend, whereas exponents 
H < 5 indicate antipersistence, i.e. past trends tend to reverse 
in the future (Fede r, 1988] >. An exponent H = | would then 
mean random uncorrelated behaviors with no stochastic mem- 
ory in time. The value H ~ 0.73 suggests that natural phenom- 
ena tend to present a persistent stochastic memory in time. 

Hurst's analysis has been applied in Astronomy mainly for 
solar phenomena. Most of the solar phenomena seem to have 
H > i, even if the actual presence of long - memory in solar 
activity is still a question of debate (Oliver & Ballester, 1998). 
Anyway, in the sun, there are also present some physical pro- 
cesses presenting different stochastic behaviors. Hanslmeier et 
al. (2000), showed that fluctuations of the velocity fields and 
the intensity fields observed in the solar upper photosphere 
have H — I . 

Let's take a practical example to illustrate Hurst's R/S anal- 
ysis. The lightcurve of a CV is a discrete record in time of 
light intensity values that usually shows both slow and rapid 
stochastic luminosity variations. Some periodic variations are 
connected to the orbital motion or to the spin of the magnetic 
white dwarf. Besides these, CVs can also present quite a vari- 
ety of irregular or quasi-periodic light oscillations. The shorter 
time scales of the flickering observed in the lightcurves of ordi- 
nary CVs are of the order of a few seconds. However, one can 
also observe smooth modulations and even flares lasting hours. 
The mean luminosity (£) T of the CV in a given time interval t 
is 

<a = -f>(0 (4) 

where t is the time coordinate. In the time lag r the accumulated 
deviation from the mean luminosity is 

X(t,T) = Y J [t(k)-(Z) T ] (5) 

k=l 

and the range R of the luminosity variations 

R(t) = [max X(t , t) - min X(t , t)] ] <,< T (6) 

As already said, the ratio R/S follows an empirical law R/S ~ 
t h like in eq 0). 

Depending on the value of H found, the behavior of the 
random fluctuations can be modeled by Wiener stochastic pro- 
cesses or, more specifically, by Fractional Brownian Motions 
(FBM). The classical Brownian motion is a particular case of 
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Fig. 1. Histogram of the Hurst exponent of all the sample of 
CVs. The distribution peaks at H ^ 0.65. Around the value 
found by Hurst for most of the natural phenomena H ~ 0.73 
there is instead a drop in the distribution. The presence of a 
second bump at H ~ 0.82 is also possible. The distribution is 
non symmetric showing a tail towards the shorter values of H. 

FBM, with H = j, without any stochastic memory in time. 
For more details, see Feder (1988). Since FBMs possess self- 
similarity, they can be also easily studied with the wavelets 
analysis (Chui 1992, Gill & Henriksen 1990, Simonsen et al., 
1998). Using the relationships between 1 /fP noises and FBMs 
( |Low en, 1999, [Gao et al., 2003| l, there is a direct relationship 
between the Hurst exponent H and the parameter a, namely 
H = 1 - or/4. 

In Fig. (fl]i is shown the histogram of the Hurst exponents 
of all the CVs analyzed by Fritz and Bruch. We see that both 
persistent and antipersistent behaviors are present. We note that 
the peak is around H ~ 0.65. At the value H ~ 0.73 found by 
Hurst in most natural phenomena there is a drastic drop. The 
presence of a second bump at H ~ 0.82 is also possible. We 
then notice that the distribution is non symmetric with a longer 
tail towards the shorter values of H. 

3. The flickering in white light of V709 Cas 

V709 Cas (RX J0028.8+5917) is a DQ Her-type system 
with orbital period of 0.2225 days ? and a 1.08M o pri- 
mary ( |Ramsay, 2000| l with a spin period P sp i„ = 312.8 s 
( |De Martino et al., 2001] >. We have analyzed six runs obtained 
in white light during a multi-site photometric campaign per- 
formed from September 22 to October 2 2000 by De Martino 
et al. (2002) with a three channel photometer placed at the 
0.8 m Tenerife (Spain), at the 1.5 m Loiano (Italy) and at the 
0.8m Beijing (China) telescopes. The integration time was in 
all cases 10s. 

The lightcurves (Fig. |2]|3]l show very rapid variations up 
to 0.2 magnitudes, indicating the presence of a strong flick- 
ering, as already reported by Kozhevnikov (2001) and De 
Martino et al. (2002). Each inset of Fig.|3]shows the "corrected 
lightcurve", i.e. the lightcurve variations with the magnitude 
rescaled with respect to zero and filtered from the fluctuations 
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Fig. 2. The full set of the original and the "corrected 
lightcurves" of V709 Cas (see text). The last night of obser- 
vations has been omitted in the right panel because of the pres- 
ence of low quality data. 



Detrended lightcurves: Power Spectrum semilog plot 
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Fig. 4. The power spectra of the six corrected lightcurves of 
V709 Cas. 
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Fig. 3. The corrected lightcurves of the six runs (a-f) of Tab. 1 . 
Time is in days and magnitudes are rescaled with respect 
to their mean magnitude. The discontinuities in the last two 
datasets are due to the presence of veils during the observa- 
tions. 



due to the orbital period of the CV. The data sets are ordered as 
in Table 1 . From up to down the first subplot (a), shows, with 
respect to the other subsets of data, a more active flickering, 
with more erratic motions, rapid and intense flares. The third 
and fifth subplot (c) and (f) show a larger variance at shorter 
timescales and shorter variance at larger timescales. They both 
present flares with larger amplitudes up to 0.2 mag. The power 
spectra of the corrected lightcurves do not show any dominat- 
ing periodicity (Fig. |4). The presence of some residual power 
(~ 10~ 4 ) is present at low frequencies, but it will not apprecia- 
bly affect the results of our analysis. 



3.1. Data analysis 

Now we determine the parameter a both with wavelets and 
Hurst R/S analysis and plot for each run the position of V709 
Cas in the parameter space a-I. with the time scale resolution of 
10 seconds. We used the C12 wavelets (Fritz & Bruch, 1998). 
To test our software we used both wavelet-generated FBMs and 
l// 8 Gaussian processes (Bak et al., 1987; Abry et al., 1995). 
Each lightcurve (a-f) was divided into series of 2" bins, with n 
is an integer varying in the interval [1,7] in order to apply the 
prescriptions of wavelet analysis. 

Tab. Q] gives, for each run, the Hurst exponent H and the 
characteristic exponent j3 of the l// 8 process obtained with 
wavelets. The table also reports the run label, the number of 
data points, the maximum and minimum of the relative mag- 
nitudes, and the percentage difference A% = |^ + ^| x 100 

between the the fractal exponent D calculated with the l// 8 
noises (D f = 1 + f) and with the R/S analysis (D H - 2- H). 

The Hurst exponents of each of the six acquisitions reported 
in Tab Q] are all larger than 1/2, indicating that the flickering 
of V709 Cas is always is a stochastic process with persistent 
memory in time. The averaged value of H is 0.63 and the stan- 
dard deviation is 0.05. The value obtained, with R/S the anal- 
ysis, reported in the table, fluctuates between H = 0.58 and 
H = 0.72. However, the differences A% are rather small. Fig. 
|5]presents the scalegrams for the six corrected lightcurves, i.e. 
the log-log plot of the variance S(t s ) vs. the chosen time scale 

ts- 

The curves rise linearly with the time intervals and remain 
quite linear for all the 6 acquisitions. This suggests that the 
flickering remains self-similar, with the same power-law ex- 
ponent, from the tens of minutes down to the 10 second time 
scale. We also see that all the curves almost overlap yielding 
a quite stable value of a, indicating that, during the 10 days- 
period of observations, the flickering and the accretion process 
of the CV remained quite stable. 
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Table 1. 



Run 


data 


Max 


min 


P 


H 


A% 


a 


1335 


0.082 


-0.075 


0.82 


0.59 


1.2 


b 


974 


0.084 


-0.093 


0.72 


0.64 


8.3 


c 


2325 


0.130 


-0.119 


0.80 


0.60 


3.8 


d 


312 


0.081 


-0.087 


0.56 


0.72 


4.3 


e 


1029 


0.108 


-0.111 


0.72 


0.64 


4.2 
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1904 
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-0.127 


0.82 


0.59 


7.3 



Scalegram of V709 Cas 
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Fig. 5. Scalegrams of the six acquisitions of the corrected V709 
Cas data. The lines corresponding to different acquisitions are 
quite similar showing a stability of a during 10 days. The small 
differences at short time scales may be due to different Poisson 
noise contributions. 

Fig- shows all the values of a and 2 of V709 Cas. The 
values span in the interval 0.9 < a < 1.8 and -1.3 < 2 < -0.6, 
typical of magnetic CVs (Bruch, 1992). In the (a, 2) plane the 
values obtained from V709 Cas show a trend. 2 is shown to 
decrease with decreasing a. This would suggest that the per- 
sistence character of the light fluctuations and the energies im- 
plied are anticorrelated. 

4. Conclusions 

We have analyzed the flickering observed in six optical photo- 
metric runs of the IP V709 Cas. In their pioneering work, Fritz 
and Bruch (1998), tried to correlate the characteristics of the 
flickering of CVs with their subtypes using wavelet analysis 
and the scalegram. The stochastic properties of flickering were 
described by using the two parameter space a-2. In this paper 
the parameter a was obtained both with the wavelets and the 
Hurst R/S analysis, the latter being a statistical robust method. 
During the 10 days of the observations the stochastic proper- 
ties of the flickering did not change indicating that the physics 
of the accretion process remained relatively stable. We also no- 
ticed that the values of a and 2 parameters, calculated at differ- 
ent epochs and at different time scales, show a trend suggesting 
that the persistence character of the light fluctuations and the 
energies implied in the accretion process are anticorrelated. In 
the (a, 2) diagram V709 Cas falls in the region of magnetic sys- 
tems. However, this CV was recognized as an IP only from its 
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Fig. 6. The a-2 parameter space for all the V709 Cas data. All 
the values for each wavelet octave are reported. The relatively 
small region of the diagram occupied by the data indicates a 
rather stable stochastic phenomenon. 



X-ray lightcurve while no spin modulation in the optical range 
is observed and the lightcurve is strongly flickering-dominated 
without any precise and stable dominating frequency, with the 
exception made for frequencies related to the orbital motion. 
Therefore, V709 Cas may represent an ideal test case where 
the a-2 classification improved by the Hurst R/S analysis can 
be used to characterize CV subtypes, especially when the clas- 
sification is uncertain. 

Hurst's R/S analysis can also provide a useful tool for un- 
derstanding the results of high speed photometry. In particular, 
it could give precious information to show down to which short 
time-scale the flickering can be still considered self-similar, 
revealing the physics behind very short time scale phenomena. 
Thus, a futuristic application might be suggested for ultra-fast 
photometry towards the quantum limit (Drav ins et al., 2 005). 
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